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Abstract 

Computations  of  divergence  and  vorticity  are  made  by  the  tri- 
angular method  for  six  100-mb  layers  (lGOO-^-00  mb)  for  those  obser- 
vation times  (ll)  when  wind  data  is  available  from  eight  stations 
in  the  BOMEX  region  during  the  period  11-28  July  1969  (Phase  IV). 
The  results  are  compared  to  cloudiness  from  satellite  pictures. 
Contemporary  correlation  coefficients  between  vorticity  and  di- 
vergence in  the  six  layers  ranged  from  -.19  to  +.26.  The  conf- 
lation in  the  lowest  layer  (1000-900  mb)  (where  Ekman  theory  pre- 
dicts a  large  (~  -1.0)  negative  correlation)  was  -.19*  The  vorticity 
and  divergence  values  had  the  same  sign  in  approximately  50^  of  the 
sample.  Despite  the  low  values  for  the  correlation  coefficients, 
the  comparison  of  the  results  to  satellite  pictures  is  encouraging 
in  that  convergence  and  cyclonic  vorticity  are  associated  with  re- 
gions of  maximum  cloudiness,  llhen   algebraic  means  are  computed 
within  each  of  the  six  layers,  the  results  indicate  the  vorticity 
has  the  same  order  of  magnitude  from  1000  to  A'-OO  mb  while  the 
divergence  decreases  by  an  order  of  magnitude  above  800  mb.  The 
divergence  and  vorticity  were  of  the  same  order  of  magnitude  from 
1000  to  800  mb.  Computation  of  mean  magnitudes  showed  no  such  de- 
crease in  the  order  of  magnitude  of  the  divergence  above  800  mb. 

All  triangles  were  subjectively  classified  as  either  disturbed, 
undisturbed  or  neither.  Comparison  of  the  mean  values  of  diver- 
gence and  vorticity  within  layers  by  the  Student's  t-test  revealed 
there  were  significant  (l$  level)  differences  between  the  vorticity 
results  for  disturbed  and  undisturbed  conditions  and  the  results 
for  divergence  indicate  significant  differences  at  the  2$  level  in 
the  lowest  layer  and  at  the  10-20;'  levels  in  the  next  two  layers 
(900-700  mb)  thereby  supporting  the  CISK  hypothesis. 
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7. 

INTRODUCTION 

The  Barbados  Oceanographic  and  Meteorological  Experiment  (BOMEX) 
was  carried  out  from  3  May  to  28  July  1969  east  of  the  island  of  Bar- 
bados (13.1  N,  59.5  W).   The  entire  experiment  was  divided  into  four 
phases  and  it  is  phase  IV  (11-28  July  1969)  with  which  this  thesis 
is  concerned.   This  phase  of  BOMEX  was  designed  to  explore  large 
convective  systems  in  the  tropical  Atlantic.   The  overall  plans, 
description  and  some  early  results  of  BOMEX  activities  are  described 
in  Davidson  (1968),  Cook  (1969),  Kuettner  and  Holland  (1969),  U.S. 
Dept.  of  Commerce  (1969),  Holland  (1970),  Friedman  and  McFadden 
(1970)  and  BOMEX  Bulletins  numbers  1-9. 

There  were  five  "stationary"  ships  plus  the  island  of  Barbados 
to  take  observations  during  phase  IV.   The  geographic  arrangement  is 
shown  in  figure  1.   In  addition,  there  were  numerous  aircraft  avail- 
able which  were  equipped  to  take  meteorological  observations.   During 
the  period  11-28  July  1969  there  were  six  synoptic-scale  disturbances 
which  passed  through  the  BOMEX  array  as  described  by  Fernandez- 
Partagas  and  Estoque  (1970),  one  of  which  was  classified  as  a  tropi- 
cal depression  (22-27  July).   The  remaining  disturbances  were  classi- 
fied as  follows:   3  weather  systems,  1  complex  weather  situation  and 
1  middle-level  weakening  cyclonic  circulation  (Fernandez-Partagas 
and  Estoque,  1970).   This  crude  method  of  classifying  tropical  dis- 
turbances suggests  a  definite  lack  of  understanding  of  the  distur- 
bance-producing mechanisms  in  the  tropics. 
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Fig.   I.       Geogrophic     orrangement     of    fixed      stot'ions       during 
BOMEX  ,  phase     TE.        (Scale  ,     1/2"  =    1°  )  . 


9. 
In  looking  for  causes  of  tropical  disturbances,  Charney  and 

Eliassen  (1964)  and  Charney  (1969)  have  proposed  conditional  insta- 
bility of  the  second  kind  (CISK)  as  a  mechanism  for  establishment 
of  low-level  water-vapor  convergence  and  a  resulting  upward  motion 
at  the  top  of  the  boundary  layer,  thus  leading  to  a  release  of  latent 
heat  above  the  boundary  layer  which  provides  the  source  of  energy 
for  maintaining  a  disturbance.   CISK  requires  an  initial  increase 
in  the  surface  relative  vorticity  as  the  initiating  impulse  to  set 
up  the  boundary-layer  convergence.   Gray  (1968)  found  a  ".  .  .strong 
association  of  tropical  disturbance  and  storm  development  with 
synoptic-scale  surface  relative  vorticity  .  .  .".   The  idea  of  CISK 
as  a  cause  of  disturbances  in  the  tropics  has  met  with  wide  approval, 
but  observational  evidence  of  such  a  mechanism  is  not  available;  in 
fact,  how  to  obtain  observations  to  confirm  the  hypothesis  is  not 
readily  apparent  to  the  author.   Sargeant  and  Ruttenberg  (1970)  have 
included  the  following  as  one  of  the  "Specific  Atmospheric  Objec- 
tives" of  an  international  tropical  experiment  which  will  probably 
be  conducted  in  1974:   "In  particular,  to  provide  sufficient  data 
to  verify  or  reject  present  formulations  of  the  CISK  mechanism,  and 
to  extend  such  formulations". 

Manabe  et,  al.  (1970)  have  attempted  to  reproduce  the  general 
circulation  features  of  the  tropics  from  a  numerical  model  for  the 
month  of  January  and  obtained  quite  good  results.   They  simulated 
the  convection  process  by  means  of  "moist-convective  adjustment", 
a  process  which  has  the  effect  of  transferring  heat  energy  from 
lower  layers  to  upper  layers,  thus  creating  a  cold  core  in  the  lower 
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part  of  the  troposphere  and  a  warm  core  in  the  upper  troposphere. 
This  process  was  initiated  in  the  model  whenever  a  layer  became 
saturated  and  the  lapse  rate  in  the  layer  became  super-moist 
adiabatic. 

In  an  attempt  to  extend  present  ideas  on  tropical  connective 
phenomena,  it  was  decided  to  use  the  BOMEX  data  to  look  at  distur- 
bances in  some  detail.   Although  the  distance  between  stations  (380- 
530  km)  precludes  investigation  of  small-scale  activity,  it  does  pro- 
vide the  necessary  data  for  investigation  of  the  disturbance-scale 
(i.e.,  sub-synoptic  scale)  as  well  as  for  comparing  differences  in 
conditions  between  disturbed  and  undisturbed  situations.   The  rela- 
tion of  the  mesoscale  structure  of  the  tropics  in  disturbed  vs.  un- 
disturbed conditions  is  now  being  recognized  as  one  of  the  most 
crucial  problems  to  solve  in  the  tropics  (Zipser,  1970). 

The  array  of  ship  and  land  stations  in  the  BOMEX  area  is  ar- 
ranged conveniently  for  application  of  the  triangular  method  of  com- 
puting divergence,  vertical  motion  and  vorticity  as  proposed  by 
Bellamy  (1949).   A  modification  of  Bellamy's  method  was  proposed  by 
Graham  (1953)  but  since  the  data  stations  in  BOMEX  are  "fixed"  the 
Bellamy  method  is  appropriate.   The  results  thus  obtained  are  most 
likely  to  be  representative  of  the  centroid  of  the  triangle  (Bellamy, 
1949).   Figure  2  shows  the  breakdown  of  the  BOMEX  area  into  eleven 
triangles  with  centroids  marked  by  <5D  .   Two  stations  have  been 
added  to  expand  the  area  of  coverage  for  which  data  is  available 
twice  daily.   The  major  disadvantage  of  the  Bellamy  method  is  the 
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Fig.   2 .      Triangular     divisions    with    expanded      BOMEX      array. 

Numbers     at    centroids     correspond    to    numbers    assigned 
to    that      triangle     for    identi  fication  .(  Scale  :   1/2"=    1°  ) 
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assumption  of  a  linear  variation  of  the  wind  between  vertices  of  the 
triangle,  an  assumption  which  may  yield  misleading  results  in  cer- 
tain cases  as  shown  by  Zipser  (1965).   Successful  application  of 
this  method  has  been  made  by  Byers  and  Rodebush  (1948)  and  Byers 
and  Hull  (1949)  in  thunderstorm  studies.   Lateef  (1967)  computed 
divergence,  vertical  motion  and  vorticity  fields  over  a  grid  cover- 
ing the  Caribbean  for  a  3-day  period  during  which  a  low-level 
easterly  wave  passed  through  the  grid.   He  used  the  pressure- 
differentiated  form  of  the  continuity  equation  to  compute  vertical 
velocities  to  offset  the  difficulty  of  obtaining  a  large  non-zero 
value  for  the  net  divergence  in  a  column.   Various  other  authors 
who  have  made  computations  of  divergence,  vertical  motion  and  vor- 
ticity in  the  tropics  utilizing  the  continuity  equation  include 
Landers  (1955) ,  Arnason  (1955) ,  Endlich  and  Mancuso  (1963) ,  Kyle 
(1970)  and  Arnason  et.  al.  (1963).   The  computations  of  Landers 
(1955)  utilized  the  modified  Bellamy  method  proposed  by  Graham 
(1953)  while  the  method  used  by  Arnason  et.  al.  (1963)  consisted 
of  extracting  wind  components  at  gridpoints  from  streamline  and 
isotach  analyses.   Kyle  (1970)  computed  vertical  motion  using  a 
latitude-longitude  grid  of  monthly  mean  values  of  the  wind  compo- 
nents.  Endlich  and  Mancuso  (1963)  performed  their  calculations  by 

the  triangular  method  on  triangles  varying  in  area  from  22000  to 

2 
920000  km  in  the  Caribbean  area.   In  contrast,  the  triangles  shown 

2 

in  figure  2  range  in  size  from  approximately  71000  to  105000  km  . 
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APPLICATION 

In  order  to  be  able  to  synoptically  analyze  the  results,  only 
those  times  when  winds  for  all  eight  stations  were  available  were 
used.   This  limited  the  data  to  00Z  and  12Z  observations  because 
that  is  the  only  data  available  for  Raizet  and  Chaguaramas  which 
is  found  in  the  Northern  Hemisphere  Data  Tabulations  (NHDT) .   This 
means  a  possibility  of  18  days  at  twice  per  day  which  would  yield 
36  data  sets.   However,  so  much  of  the  ship  wind  observations 
were  missing  that  the  end  result  was  11  data  sets,  and  they  aren't 
all  quite  complete.   The  observation  times  finally  selected  for 
this  study  are:  14/OOZ,  14/12Z,  18/00Z,  20/00Z,  21/12Z,  22/12Z, 
25/12Z,  26/00Z,  26/12Z,  28/OOZ  and  28/12Z.   Tables  3  to  5  in  Appen- 
dix B  summarize  the  identification  of  the  triangles  and  the  loca- 
tions of  the  centroids. 

In  order  to  interpret  the  results  in  comparison  with  disturbed 
or  undisturbed  situations,  a  method  of  classifying  the  triangles 
as  disturbed  or  undisturbed  is  needed.   Such  a  method  is  available 
by  means  of  enlargements  of  ATS-III  and  ESSA-9  satellite  pictures. 
By  tracing  the  cloud  cover  in  the  BOMEX  region  onto  acetate  and  then 
projecting  the  result  with  an  overhead  projector,  enlargements  of 
the  satellite  picture  can  be  obtained  to  almost  any  scale.   For  this 
study  the  pictures  were  enlarged  to  a  scale  of  1"  =  1°  of  latitude 
(results  presented  later  are  at  a  scale  of  1/2"  =  1°  of  latitude). 
With  a  copy  of  figure  2  at  the  same  scale  superimposed  on  the  copy 
of  the  satellite  picture,  a  subjective  classification  of  U  (undis- 
turbed) ,  D  (disturbed)  or  N  (neither  disturbed  nor  undisturbed)  was 
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made  for  each  triangle  at  each  observation  time  depending  on  the  areal 
extent  of  cloud  cover  for  each  triangle.   Allowances  were  made  for 
cases  when  the  satellite  picture  was  not  at  the  same  time  as  the  ob- 
servations.  The  criteria  used  for  the  classifications  were  as 
follows:  U,  less  than  10%  of  area  covered  by  clouds;  N,  less  than 
50%  of  area  covered  by  clouds  but  more  than  10%;  D,  50%  or  more  of 
area  covered  by  clouds.   Out  of  a  total  of  121  such  classifications 
(11  triangles  for  each  of  the  11  observation  times  used)  there  were 
52  =  U,  27  =  N  and  42  =  D.   The  results  of  this  classification  are 
listed  in  Appendix  B,  table  6. 

Appendix  A  describes  the  method  of  computing  divergence  and 
vorticity.   In  addition  to  the  triangle  parameters,  the  data  needed 
for  the  computations  are  the  wind  direction  and  wind  speed  profiles 
for  each  of  the  stations.   For  this  study  the  winds  were  obtained 
at  50-mb  intervals  from  1000  to  400  Bib  and  then  vectorially  averaged 
by  100-mb  layers  according  to  the  trapezoidal  rule,  thus  resulting 
in  six  layers,  each  containing  a  vector-mean  wind.   The  vector- 
mean  winds  were  then  used  to  compute  the  mean  divergence  and  mean 
vorticity  for  each  of  the  six  layers  of  each  of  the  eleven  triangles. 
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DATA  UTILIZED 

The  data  utilized  for  the  computations  were  the  u  and  v  wind 
components  at  50-mb  intervals  from  1000  to  400  mb.   The  data  for 
Barbados,  Raizet  and  Chaguaramas  were  used  exactly  as  given  in 
the  NHDT.   Data  for  the  five  "stationary"  ships  were  complicated 
by  the  fact  that  the  ships  were  actually  not  stationary  due  to  the 
failure  of  the  mooring  lines  for  all  five  ships  before  the  begin- 
ning of  phase  IV  (BOMAP,  1971).   Since  it  was  planned  to  use  the 
processed  rawinsonde  observations  from  the  "A0"  (BOMEX  Bulletin 
No.  9,  1971)  data  reduction  process  performed  at  NASA's  Mississippi 
Test  Facility,  this  presents  the  problem  of  correcting  the  winds 
for  ship's  motion.   A  search  through  ship's  logs  and  course  and 
speed  tables  provided  by  BOMAP  (BOMAP,  1971)  revealed  that  there 
were  several  inconsistencies  in  regards  to  the  necessary  correction, 
(For  example,  ship's  course  before,  during  and  after  could  not  be 
verified  when  compared  to  gyro  headings  found  in  the  BOOM  data 
[BOMAP,  1971]).   Fortunately,  however,  it  appears  that  this  cor- 
rection is  very  small  (because  tabulated  ship  speeds  are  small)  for 
all  ships  except  ROCKAWAY.   The  tabulated  ship  speeds  for  ROCKAWAY 
were  as  large  as  10. 5  knots.   Due  to  the  uncertainty  in  such  a 
large  correction  it  was  decided  to  use  the  data  sent  by  ROCKAWAY 
via  teletype  with  the  idea  that  the  observers  on  ROCKAWAY  made  the 
necessary  corrections  to  the  data  before  sending  it. 

An  additional  problem  with  ROCKAWAY  was  the  inability  to  ac- 
quire the  balloon  for  6  to  10  minutes  after  launch  (BOMAP,  1971), 
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as  evidenced  by  the  missing  winds  in  the  lowest  5000  feet  or  so. 
To  alleviate  this  problem  the  teletype  winds  were  broken  down  into 
u,  v  components  and  plotted  on  a  vertical  logarithmic  scale  using  a 
standard  tropical  atmosphere  (Schacht,  1946;  and  Riehl,  1954).   The 
winds  at  the  lowest  reported  level  were  then  linearly  connected  to 
the  surface  wind  and  the  necessary  winds  were  read  off  at  50-mb 
intervals. 

Except  for  specific  cases  to  be  noted  later,  the  "A0"  data 
was  used  for  the  remaining  ships  (at  present  there  is  no  choice 
but  to  use  it  for  Mount  Mitchell).   The  u,  v  winds  were  taken  di- 
rectly from  the  "A©11  data  at  50-mb  intervals.   In  a  few  cases 
winds  were  missing  at  some  of  the  50-mb  intervals  (in  particular, 
for  DISCOVERER)  and  in  such  cases  the  needed  data  was  obtained  by 
applying  techniques  similar  to  that  described  above  for  ROCKAWAY. 
The  only  corrections  made  to  this  data  were  for  RAINIER  which  had 
a  voltage  reference  problem.   Comment  cards  acquired  from  BOMAP 
on  the  performance  of  each  individual  sounding  indicate  that  the 
computed  wind  directions  for  RAINIER  are  in  error  by  140-220°.   This 
is  particularly  evident  if  wind  directions  computed  from  "Ao"  data 
are  compared  to  wind  directions  reported  in  the  teletype  data.  Con- 
sequently the  error  reported  by  BOMAP  was  applied  to  all  the  wind 
directions  for  RAINIER  as  computed  from  the  "A0"  data.   Even 
though  no  correction  was  tabulated  for  RAINIER  at  26/12Z,  a  cor- 
rection of  140°  was  applied  based  on  soundings  at  times  before  and 
after  26/12Z  and  based  on  comparisons  with  the  teletype  data. 
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Teletype  data  was  also  substituted  for  "A0"  data  in  the  following 

three  cases  for  the  reasons  specified: 

21/12Z  Oceanographer  ;  Many  winds  were  missing  in  the  "A0" 
data  but  complete  wind  data  was  available  in  teletype  form. 

22/12Z  Oceanographer;  There  were  no  winds  in  the  "A0"  data 
but  complete  wind  data  was  available  in  teletype  form. 

22/12Z  Discoverer;  The  "A0"  winds  were  missing  at  many 
levels  but  complete  wind  data  was  available  in  teletype  form. 
There  was  no  data  at  all  for  OCEANOGRAPHER  at  28/12Z  but  there  was 
MA0"  data  for  28/10Z  and  28/1343Z.   In  order  to  complete  the  data 
set  at  28/12Z,  the  winds  at  the  two  off  times  were  averaged  and 
used  in  the  computations  to  obtain  an  approximation  to  the  values 
for  28/12Z. 

A  final  comment  regarding  corrections  to  the  island  data  is 
necessary.   No  terrain  effects  have  been  considered  for  Barbados, 
Raizet  and  Chaguaramas.   As  pointed  out  by  Lilly  and  LaSeur  (1956) 
the  nature  of  a  correction  at  Chaguaramas  is  very  complicated. 
This  was  also  discussed  by  Zipser  (1965)  who  made  "ad  hoc"  correc- 
tions to  the  Chaguaramas  winds.   Zipser' s  correction  to  the  winds 
at  Raizet  amounted  to  increasing  the  reported  mean  wind  speed 
(surface  -  850  mb)  by  17%.   For  this  study,  no  corrections  have 
been  applied  to  the  reported  winds  for  any  of  the  island  stations. 
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RESULTS 

Some  of  the  results  of  this  study  are  presented  in  this  section 
in  the  form  of  synoptically  analyzed  charts  of  the  BOMEX  region. 
Vorticity  and  divergence  for  layers  1  and  2  (1000-900  mb  and  900- 
800  mb)  for  five  synoptic  times  are  presented  in  the  figures  3-12. 
Results  for  all  synoptic  times  used  and  all  layers  are  presented 
in  tables  7-12  in  Appendix  B.   In  addition,  figures  3-12  show  de- 
pictions of  cloud  cover  obtained  from  ATS-III  and  ESSA-9  satellite 
pictures  for  comparison  to  the  analyzed  results  which  are  presented. 
These  depictions  were  obtained  from  the  satellite  pictures  as  des- 
cribed in  the  Application  section. 

The  figures  presented  in  this  section  cover  two  specific 
cases:  1)   14  July  1969;  Weather  system  (Fernandez-Partagas  and 
Estoque,  1970)  which  moved  through  the  northern  half  of  the  BOMEX 
region  from  east  to  west  and  2)  25-26  July  1969;  Tropical  depression 
(Fernandez-Partagas  and  Estoque,  1970)  which  emerged  from  the  ITCZ 
and  moved  northwest  through  the  BOMEX  region. 

A  comparison  of  the  satellite  picture  depiction  at  14/00Z  (keep 
in  mind  the  satellite  picture  is  for  13/1817Z)  with  the  analysis  of 
divergence  and  vorticity  shows  nothing  very  striking,  perhaps  due 
to  the  relative  inactivity  in  the  BOMEX  region.   Perhaps  the  most 
noticeable  feature  is  the  convergence  which  seems  to  be  coming  into 
the  region  from  the  east.   At  14/12Z  the  convergence  has  extended 
slightly  farther  into  the  BOMEX  region  and  appears  to  be  somewhat 
associated  with  the  region  of  cloudiness  which  is  approaching  from 
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Fig.     3.    Anolyzed     results    of     divergence     and 
vorticity      for      I4/00H  ,     layer      I     ,    with 
cloud      cover        for     13/1817  2- 
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Fig.    4'.      Same     os     figure       3     except    layer     2. 
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Fig.  5:  Analyzed  results  of  divergence  and  vortici'y 
for  14/12  2  ,  layer  I  ,  with  cloud  cover 
for     14/  11362. 
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Fig.  6  '.     Same    as    figure    5    except     layer      2 
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the  east.   Likewise  there  is  a  region  of  cyclonic  vorticity  which 
appears  to  be  associated  with  the  cloud  system  and  anticyclonic  vor- 
ticity out  ahead  of  it. 

Comparing  the  25/12Z  results  (the  satellite  picture  reproduc- 
tion is  for  25/1532Z)  shows  the  ITCZ  to  be  associated  with  the  region 
of  strongest  convergence.   The  stipled  regions  in  the  cloud  depic- 
tion correspond  to  regions  on  the  picture  which  were  very  bright. 
Comparison  with  the  vorticity  analysis  shows  the  ITCZ  region  to  be 
associated  with  the  cyclonic  vorticity;  however,  the  orientation  of 
the  vorticity  lines  appears  to  be  transverse  to  the  ITCZ. 

For  comparison  of  results  at  26/00Z  it  becomes  necessary  to  use 
a  25/1826Z  satellite  picture.   Unfortunately  there  isn't  much  dif- 
ference between  the  25/15Z  and  25/18Z  pictures  but  the  difference 
in  the  results  for  divergence  and  vorticity  is  considerable.   The 
region  of  strongest  convergence  is  still  associated  with  the  region 
of  maximum  cloudiness  but  now  there  is  also  a  strong  (^4  x  10   ) 
center  of  cyclonic  vorticity  which  appears  to  be  located  in  the 
ITCZ  region. 

At  26/12Z  the  pattern  has  changed  considerably.   There  is  a 
broad  expanse  of  cloud  cover  which  covers  almost  the  entire  BOMEX 
region.   There  is  a  strong  (-3.5  x  10   )  center  of  convergence  in 
the  center  of  the  region  and  the  small  region  of  divergence  (13N, 
61W)  agrees  well  with  a  clear  region  on  the  satellite  picture.   The 
vorticity  results  are  not  as  good  except  for  the  cyclonic  vorticity 
region  in  the  northern  half  of  the  BOMEX  region. 
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Fig.    7:    Analyzed     results    of    divergence     and    vorticity 
for     25/122    ,    layer    I    ,    with     cloud     cover 
for      25/15322- 
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Fig.  8'.     Same     as    figure      7     except      layer     2 
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Fig.  9*    Analyzed    results     of    divergence    ond     vorticity 
for     26/00  2  ,     layer    I  ,    with    cloud    cover 

for     25/I826Z. 
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Fig.    10 :      Same     as    figure      9    except    layer      2 
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Fig.    II  '     Analyzed     results   of  divergence    and    vorticity 
,     for    26/122  ,  layer      I,  with     cloud    cover 
for     26/1233?. 
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Fig.   12"      Same    as     figure      II     except     layer      2 
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In  order  to  compare  these  results  with  those  predicted  by  Ekman 
theory,  a  correlation  coefficient  between  divergence  and  vorticity 
was  computed.   Correlations  were  computed  for  all  layers  to  look 
at  differences  between  layers.   To  minimize  the  effects  of  round- 
off error,  all  values  were  multiplied  by  10  before  computing  any 
means  or  sums.   The  results  of  this  are  as  follows  (%  with  same 
sign  means  percent  of  sample  where  vorticity  and  divergence  had 
same  sign) : 

Layer  (mb)     Sample  Size    %  Same  Sign     Correlation  Coef- 
ficient 


1000-900 

117 

43.6 

-.19 

900-800 

117 

55.5 

-.03 

800-700 

113 

58.4 

+  .26 

700-600 

113 

56.6 

+  .15 

600-500 

109 

49.5 

-.05 

500-400 

96 

50.0 

+  .08 

These  values  in  the  lower  layers  are  somewhat  disappointing  from 
what  one  would  expect  from  Ekman  theory.   Nevertheless  it  must  be  re- 
membered that  this  is  a  linear  correlation  coefficient  and  as  such  is 
merely  a  measure  of  the  scatter  about  a  best-fit  line  of  the  values  of 
divergence  and  vorticity.   Another  possibility  is  that  there  is  a  lot 
of  "noise"  in  the  data  which  was  used.   It  was  anticipated  that  by 
using  layer-mean  winds  in  the  computations,  instead  of  winds  at  dis- 
crete levels,  much  of  the  unwanted  noise  or  inconsistencies  would  be 
minimized.   Perhaps  the  mixing  of  "A0"  data  with  teletype  data  as  well  as 
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linearly  adjusting  the  ROCKAWAY's  low-level  winds  contributed  to 
data  noise.   The  low  values  for  the  correlation  coefficient  could 
have  been  anticipated  from  the  higher-than-expected  percentage 
(50%)  of  cases  where  the  divergence  and  vorticity  had  the  same  sign. 

The  possibility  of  computing  divergence  and  vorticity  on  one 
scale  with  data  contaminated  by  local  small-scale  winds  may  also  be 
a  problem.   It  merely  means  the  assumption  of  a  linear  change  in  the 
wind  between  vertices  of  the  triangle  does  not  apply.   By  using  layer- 
mean  winds,  however,  the  effects  of  such  wind  regimes  should  be 
minimized. 

The  profiles  of  the  mean  values  of  divergence  and  vorticity 
provide  an  interesting  result  which  is  summarized  in  table  1.   (The 
six  samples  which  were  used  to  obtain  the  results  in  table  1  are 
contained  in  Appendix  B,  tables  7-12).   The  significant  fact  which 
is  immediately  obvious  is  that  the  mean  (algebraic)  divergence  is 
an  order  of  magnitude  greater  below  800  mb  than  it  is  above  while 
the  magnitude  of  the  mean  (algebraic)  vorticity  is  about  the  same 
from  1000  to  400  mb.   The  sample  size  is  not  very  large  even  in 
the  lowest  layers  where  there  is  more  data  available,  but  the  maxi- 
mum values  of  cyclonic  vorticity  and  convergence  do  appear  to  occur 
in  the  lowest  layers.   Furthermore,  in  the  mean,  the  vorticity  and 
divergence  are  of  opposite  sign  and  are  of  the  same  order  of  mag- 
nitude from  1000  to  800  mb . 

Table  1  also  contains  the  values  of  the  mean  magnitude  of 
divergence  and  vorticity  within  each  of  the  six  layers.   When  looked 
at  in  this  manner  the  magnitude  of  the  divergence  does  not  decrease 
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by  a  factor  of  10  above  800  mb.  Comparison  of  these  results  with 
those  of  Endlich  and  Mancuso  (1963,  p.  21)  show  a  tendency  for  the 
magnitude  of  the  vorticity  to  be  5-20$  larger  than  theirs  up  to 
500  mb  while  the  magnitude  of  the  divergence  is  20-60$  larger  than 
theirs.  Above  500  mb  the  magnitudes  of  both  divergence  and  vorticity 
given  in  table  1  are  less  than  those  of  Endlich  and  Mancuso  (1963). 

Similar  computations  of  algebraic  means  and  mean  magnitudes 
were  carried  out  for  each  of  the  three  categories  of  degree  of 
disturbance  (U,  N  and  D)  and  the  results  are  presented  in  table  2. 
Although  this  results  in  three  small  samples,  particularly  for 
category  N,  the  results  between  categories  do  appear  to  differ 
dramatically  in  several  respects.  The  mean  (algebraic)  vorticity 
is  quite  large  and  cyclonic  from  1000  to  600  mb  in  category  D 
whereas  it  is  weak  and  anticyclonic  in  category  U.  The  mean 
vorticity  is  also  cyclonic  from  1000  to  600  mb  in  category  N  but 
not  as  large  as  in  category  D.  The  variance  of  the  vorticity  in 
category  D  is  quite  large  and  increases  from  1000  to  700  mb  even 
though  the  algebraic  mean  remains  practically  constant.  The  mean 
(algebraic)  divergence  also  differs  considerably  between  categories 
U  and  D.  In  category  U  there  is  very  slight  convergence  but  it  is 
an  order  of  magnitude  less  than  the  convergence  from  1000  to  800  mb 
in  category  D.  The  results  for  category  N  appear  to  be  intermediate 
to  those  of  categories  U  and  D  (as  would  be  expected)  except  for  the 
convergence  in  layer  3* 

To  determine  if  the  results  in  table  2  show  any  significant 
differences  between  categories,  the  Students  t-test  (Panofsky  and 
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Category 

U 

Layer 
(mb) 

Sample 
Size 

? 

Vari- 
ance 

\IL 

W 

Vari- 
ance 

|VAV| 

1000-900 

51 

-.05 

1.12 

.89 

-.08 

.93 

.74 

900-800 

51 

-.26 

1.76 

1.11 

0.00 

1.09 

.83 

800-700 

50 

-.13 

1.09 

.83 

-.04 

1.12 

.83 

700-600 

50 

.06 

1.20 

.89 

.01 

1.72 

1.00 

600-500 

49 

-.19 

1.37 

.87 

.08 

2.01 

1.01 

500-400 

43 

-.75 

1.10 

1.01 

-.13 

1.19 

.89 

Category  N 


1000-900 

24 

.56 

1.06 

.98 

-.31 

1.67 

1.04 

900-800 

24 

.26 

1.43 

.92 

-.24 

2.92 

1.33 

800-700 

23 

.12 

1.70 

.90 

-.26 

2.39 

1.17 

700-600 

23 

.12 

2.44 

1.01 

.02 

1.53 

.82 

600-500 

23 

-.23 

1.3^ 

.93 

.38 

1.45 

.92 

500-400 

21 

-.20 

.90 
Category 

.84 
D 

.13 

1.01 

.83 

1000-900 

42 

.94 

2.44 

1.46 

-.61 

1.43 

1.05 

900-800 

42 

.98 

2.51 

1.5^ 

-.27 

1.93 

1.08 

800-700 

40 

.98 

3.10 

1.66 

.17 

.96 

.78 

700-600 

40 

.66 

2.69 

1.41 

0.00 

1.10 

.84 

600-500 

37 

-.02 

1.92 

1.13 

-.08 

1.36 

.93 

500-400 

32 

-.20 

1.01 

.75 

-.04 

.94 

.71 

Table  2.  Summary  of  mean  values  of  divergence  and  vorticity     ,5 
within  categories  of  degree  of  disturbance.  Units  are  10  , 
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Brier,  1958)  was  performed  on  the  computed  means  of  vorticity  and 
divergence  by  layers  between  the  three  categories.  If  the  results 
are  significant  they  would  be  expected  to  occur  between  categories 
U  and  D. 

The  results  of  this  test  were  as  follows: 
Categories  U  vs.  D 


Layer 

^ 

$  Chance 
that  U=D 

~£y»</ 

$  Chance 
that  U=D 

1 

3-574 

<1 

2.333 

2 

2 

4.067 

<1 

1.071 

10-20 

3 

3.662 

<1 

.932 

10-20 

k 

2.073 

<3 

.024 

>45 

5 

.617 

25-30 

.538 

25-30 

6 

2.204 

<3 
Categorie 

.363 
c  U  vs.   N 

30-40 

Layer 

k 

$  Chance 
that  U=N 

~fcy'W 

i>  Chance 
that  U=M 

1 

2.313 

O 

.837 

20-25 

2 

1.591 

5-10 

.753 

20-25 

3 

.859 

10-20 

.702 

20-25 

4 

.194 

40-45 

.023 

>45 

5 

.124 

45 

.884 

10-20 

6 

1.912 

2-5 

.918 

10-20 
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Categories  N  vs.  D 


Layer 

1.041 

$  Chance 
that  N=D 

10-20 

.941 

io  Chance 
that  N=D 

1 

10-20 

2 

1.926 

2-5 

.074 

>45 

3 

2.005 

1-2 

1.313 

5-10 

4 

1.270 

10-20 

.046 

>45 

5 

.596 

25-30 

1.^51 

5-10 

6 

.029 

>45 

.619 

25-30 

The  results  are  evidently  very  significant.  During  disturbed 
conditions  there  is  significantly  more  cyclonic  vorticity  than 
during  undisturbed  conditions.  Furthermore  there  is  significantly 
more  convergence  in  the  lowest  layer  in  disturbed  versus  undisturbed 
conditions  and  the  results  for  convergence  in  the  900  to  700-mb 
layers  could  easily  be  called  meteorologically  significant.  The 
results  between  categories  U  and  N  and  categories  N  and  D  occassion- 
ally  show  significant  differences  but  this  is. probably  caused  by  the 
N-sample  being  a  mixture  of  both  U  and  D  values  as  well  as  being 
a  small  sample.  It  is  interesting  to  note  that  there  is  apparently 
no  significant  difference  between  samples  for  convergence  in  the 
700  to  600-mb  layer  because  all  three  samples  were  essentially 
non-divergent  according  to  the  algebraic  mean  (see  table  2). 

The  results  between  categories  U  and  D  are  in  direct  support 
of  the  CISK  hypothesis  (Charney  and  Eliassen,  1964)  in  that  there 
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is  increased  cyclonic  vorticity  and  increased  convergence  in  the 
lowest  layers  during  disturbed  conditions.  In  addition  the  results 
presented  here  indicate  the  increased  cyclonic  vorticity  is 
present  up  to  600  mb  whereas  the  increased  convergence  occurs  up 
to  800  mb. 
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CONCLUSIONS 

From  the  results  presented  in  this  paper,  it  appears  that  the 
triangular  method  of  computing  divergence  and  vorticity  in  the 
tropics  yields  acceptable  results.  The  contemporary  correlations 
of  vorticity  and  divergence  in  the  lowest  layers  were  disappointing, 
however  the  good  agreement  which  was  obtained  when  comparing  sate- 
llite pictures  to  analyzed  results  was  encouraging.  When  the  results 
were  divided  into  categories  depending  on  whether  they  occurred  dur- 
ing disturbed  or  undisturbed  conditions  as  determined  by  cloud  cover 
from  satellite  pictures,  there  were  significant  differences  between 
disturbed  and  undisturbed  conditions  for  both  vorticity  and  diver- 
gence. There  was  a  significant  increase  in  cyclonic  vorticity  and 
convergence  which  is  in  support  of  the  CISK  hypothesis  (Charney  and 
Eliassen,  1964) .  Algebraic  mean  values  of  divergence  and  vorticity 
for  the  entire  sample  had  the  same  order  of  magnitude  below  800  mb 
but  the  divergence  decreased  by  an  order  of  magnitude  above  800  mb. 
Furthermore  the  mean  (algebraic)  vorticity  retained  the  same  order  of 
magnitude  from  1000  to  400  mb. 

It  is  suggested  that  additional  computations  be  made  of  the 
nature  described  in  the  previous  sections  but  without  the  emphasis 
on  those  times  when  all  eight  stations  reported.  Rather  the  empha- 
sis should  be  shifted  so  as  to  obtain  a  very  large  sample  of  data 
to  determine  if  the  results  in  tables  i  and  2  remain  valid.  A  more 
objective  method  of  classifying  the  results  into  categories  of  undis- 
turbed or  disturbed  is  needed.   In  addition  it  is  not  possible  from  this 
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study  to  conclude  the  reason  for  the  divergence  and  vorticity  having 
the  same  sign  50%  of  the  time  in  the  lowest  layers,  a  result  which 
is  contradictory  to  Ekman  theory. 

Classification  of  the  triangular  regions  as  disturbed  or  un- 
disturbed might  be  improved  by  combining  the  method  used  in  this 
study  with  a  method  which  is  more  objective.   Such  a  method  is  avail- 
able using   ©^  >  the  equivalent-potential  temperature  (Rossby,  1932 
and  Hess,  1959).   Garstang  et.  al.  (1967)  found  that  higher  values 
of   C7g  occurred  at  all  levels  when  the  atmosphere  was  disturbed. 
In  some  preliminary  reduction  of  the  BOMEX  data  for  12-15  and 
23-28  July  1969  (  168  soundings)  the  author  found  the  vertically- 
averaged  value  of  @e  (surf  ace  to  400  mb)  to  vary  from  325.  OK  to 
342. IK,  the  higher  values  occurring  during  disturbed  conditions. 
The  advantage  of  using  ^~  is  that  it  combines  the  effects  of  tem- 
perature and  moisture  into  one  parameter.   Furthermore  such  an  ob- 
jective method  could  be  used  by  itself  for  times  when  no  satellite 
pictures  are  available.   Such  an  approach  might  yield  results  which 
show  more  significant  differences  between  the  mean  divergence  and/or 
vorticity  in  disturbed  and  undisturbed  conditions.   Significant 
differences  (5%  level  according  to  Student's  t-test)  were  not 
found  in  this  study  although  the  differences  between  the  undis- 
turbed versus  disturbed  divergence  in  the  lowest  three  layers  were 
significant  at  15-35%  levels. 

Due  to  the  difficulties  encountered  with  the  "Ao"  data  format, 
it  would  seem  that  the  teletype  data  is  still  the  best  (most  reliable) 
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data  source,  even  though  it's  been  almost  two  years  since  the  ter- 
mination of  BOMEX.   Since  the  observers  on  board  the  ships  are  "pro- 
fessionals" at  their  jobs,  it  seems  reasonable  that  any  corrections 
which  need  to  be  made  to  the  data  (e.g.  ship  motion)  would  be  made 
by  the  observers  before  sending  the  results  of  a  launch.   It  is 
unfortunate  that  detailed  records  are  not  readily  available  which 
reflect  what  corrections,  if  any,  were  made  by  the  observers  be- 
fore sending  the  teletype  message. 

Since  no  corrections  were  made  to  the  winds  at  the  island  sta- 
tions, it  is  suggested  that  a  method  for  adjusting  the  winds  repor- 
ted at  these  stations  be  adopted  before  conducting  additional  re- 
search with  triangles  affected  by  these  stations.   Corrections  were 
suggested  by  Zipser  (1965)  for  Raizet  and  Chaguaramas; however  the 
"correction"  used  by  him  at  Chaguaramas  was  "ad  hoc". 
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APPENDIX  A 
The  Computation  Procedure 

Numerous  triangles  can  be  formed  within  the  BOMEX  array  as 
shown  in  figure  2  which  can  be  used  to  calculate  horizontal  con- 
vergence (  Vh*^  )  as  shown  by  Bellamy  (1949): 


wh 


ere   J\/;|   is  the  observed  wind  speed,  h.  =  the  altitude  of  the 


triangle  through  vertex  i,  d.  =  the  observed  wind  direction  and 


A. 


■^iK  =  the  azimuth  of  the  side  of  the  triangle  opposite  vertex  i. 
With  ]V;(  in  (m  sec   )  and  h.  in  (m)  ,  the  units  of  7n*V  become  sec 
Bellamy  (1949)  describes  construction  of  an  overlay  to  make  the 
necessary  computations  or,  as  an  alternative,  to  construct  tables 
which  can  be  entered  with  wind  direction  aid  wind  speed.   However, 
with  the  advent  of  the  electronic  computer  direct  computations  using 
(1)  are  simplified. 

Computations  of  vorticity  can  also  be  performed  using  (1)  if 
d.  is  increased  by  90°  (Bellamy,  1949),  such  that 


J  ;  - 1    U  • 


c--«   Kl  (2) 

Since  data  were  used  at  constant  pressure  surfaces,  the  compu- 
tation of  vertical  motion  with  the  continuity  equation  is  straight- 
forward. 

%££  -  _v\V  (3) 
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Integrating   leads   to 

oo^co,  -  7-VAp  (A) 

/</vvvv 

where  the  symbol  v    )  implies  a  vertical  average. 

Using  the  boundary  condition  that 

CO  =  0   at   p  -  1000  mb  (5) 

enables  computation  of  vertical  velocities  by  successive  integra- 
tions of  (3) . (  A  review  of  various  methods  of  computing  vertical 
motion  are  contained  in  Panofsky  [1946]).   An  additional  assump- 

A'VVW*  A'-'.'i 

tion  that  V*^V  ~  ^  *'V    was  also  made. 
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APPENDIX  B 
Table  3.   Latitude  and  longitude  of  data  stations: 

Station  Latitude  (N)      Longitude  (W) 


1. 

Oceanographer 

7.5 

52.7 

2. 

Mount  Mitchell 

10.5 

56.5 

3. 

Discoverer 

13.0 

54.0 

4. 

Rockaway 

15.3 

56.6 

5. 

Rainier 

17.5 

54.0 

6. 

Barbados 

13.1 

59.5 

7. 

Raizet 

16.3 

61.5 

8. 

Chaguaramas 

10.7 

61.6 

Table  4.   Designation  of  triangles: 

Triangle  No.  Vertices  of  the  triangle 

1  Oceanographer,  Chaguaramas,  Mount  Mitchell 

2  Oceanographer,  Mount  Mitchell,  Discoverer 

3  Mount  Mitchell,  Chaguaramas,  Barbados 

4  -  Mount  Mitchell,  Barbados,  Discoverer 

5  Discoverer,  Barbados,  Rockaway 

6  Mount  Mitchell,  Barbados,  Rockaway 

7  Mount  Mitchell,  Rockaway,  Discoverer 

8  Discoverer,  Rockaway,  Rainier 

9  Rockaway,  Raizet,  Rainier 

10  Barbados,  Raizet,  Rockaway 

11  Chaguaramas,  Raizet,  Barbados 
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APPENDIX  B 
Table  5.   Latitude  and  longitude  of  triangle  centroids: 


Triangle  No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 


Latitude(N) 

Longitude (W) 

9.5 

56.9 

10.3 

54.4 

11.4 

59.2 

12.2 

56.6 

13.8 

56.7 

13.0 

57.5 

12.9 

55.7 

15.2 

54.9 

16.3 

57.4 

14.9 

59.2 

13.4 

60.9 
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